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TECHNICAL MEMORANDUM 


COMPONENT RESPONSE TO RANDOM VIBRATORY MOTION 
OF THE CARRIER VEHICLE 

SECTION 1. INTRODUCTION 


In this treatment of component response to local random vibratory motion of the 
carrier vehicle, the component plus supporting structure is modeled as the system 
shown in Figure 1. The component model is allowed two degrees-of- freedom, one 
translational and one rotational, and is excited by a random translatory motion of the 
base whose acceleration power spectral density (PSD), herein denoted by G..(f), is 

presumed known. Prescription of the base acceleration PSD is done in the manner 
indicated by the inset in Figure 2 which admits the analytical representation appear- 
ing beneath the diagram of Figure 3. 

Since the base motion is prescribed only to the extent that its acceleration PSD 
is given, the ’’time response,” i.e., a time history of the system configuration coor- 
dinates and their first and second time derivatives, is out of the question. The word 
"response” is here to be interpreted as "mean square response," that implying the 
mean squares of the system coordinates and their time derivatives pertinent to the 
frequency interval over which G«(f) is specified. 


SECTION 2. FUNDAMENTAL RELATIONS 


Whether interest lies in "time response" or "mean square response," the source 
of certain fundamental relations, necessary to computation, is the system of differential 
equations descriptive of the motion. Treating the component model as a perfectly 
rigid body, invoking Newton and the principle of angular momentum, and making the 
usual small angle approximations, the equations of motion may be written as equations 
(1) and (2). 


mx = -mg - Kj (x - 6 g^, ^ - L^ 6 - u) - (x - D^ 0 - u) 


K 2 (x - 6 ST 2 + L 2 0 “ “ C 2 ^ + D 2 ® " u) 


( 1 ) 


16 = (x - 6 Qrr 1 - L 1 0 - u) + C 1 D 1 (x - D 1 0 - u) 


ST ,1 “1 


1"1 


- K 2 L 2 ( x " 5 ct> 9 + L <> 0 " u ) “ C 9 D o (x + D 0 0 - u) 


ST, 2 2 


2 2 


( 2 ) 


Recognizing the simplifications possible via the relations (3), the conditions for static 
equilibrium , 
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one can write the matrix equivalent of equations (1) and (2) as 
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If the (y , 0) -description of system configuration is preferred to the (x, 0)- description, 
then one has only to make the substitution x = y+u in equation (4) to get 
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C 2 D 2 " C 1 D 1 
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_ C 2 D 2 - C 1 D 1 

C 1 D 1 2 + C 2 D 2 2 . 
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-mu" 

_ K 2 L 2 - K ! L 1 

K 1 L 1 2 + K 2 L 2 2 
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(5) 


The transfer functions K = x,y,0, essential to computation of the mean 

squares of the system output variables, may be found by applying the Laplace trans- 
formation to equations (4) and (5), assuming zero initial conditions, and solving for 
the transform ratios <2f(0/<^(u), £ = x,y,9. (Obviously, it is not necessary to 
apply the transformation to both equations (4) and (5) since one can choose to work 
with either equation (4) or (5), then, having found either or (y) A5f(u) , 

find the other via the relation x = y+u.) As should be expected, the expression for 
T 0 / u ( s ) as determined by equation (4) is equivalent to that determined by equation 

(5), Thus, 



( 6 ) 


T x / U (s) = { a 2 2 (s) t a 1 i(s> - ms 2 ] - a 2 2 (s) } / A(s) 

T y/u (s) = - ms 2 a 22 (s) / A(s) (7) 

T 0/U (s) = ms 2 a 12 (s) / A(s) (8) 

where 

aii(s) = ms 2 + ( Cl + C 2 ) s + K 1 + K 2 
a l2 (s > = (C 2 D 2 ' C 1 D 1> s + K 2 L 2 ' K 1 L 1 

(9) 

a 22 (s) = Is 2 + (C^ 2 + C 2 D 2 2 ) s + K 1 L 1 2 + K^ 2 

2 

A(s) = a 1;l (s) a 22 (s) - a 12 (s) 

Having found T^ u (s), £ = x,y , 9, it is an easy matter to find an expression for the 
PSD of £, G^(f), by appealing to the well known general relation (valid for linear 
systems and accepted here without dispute) 


G^Cf) = |T c/n (2 tt jf) | 2 G n (f) , (10) 

the symbol t, denoting an output quantity of a system with input n. Obviously, 

G £ (f) = l T £/ u (2 u if) I 2 G u (f) > ^ = x »y> e • (id 


But, since it is G~(f) that is prescribed, not G u (f), equation (11) will not completely 
define G^(f) until an expression for G u (f) is found. To that end one can write 


u/ii 


..(s) = <^(u)/^(u) = 


S 2 ^(u) 
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which , in conjunction with equation ( 10) , yields 


G u ( f ) 


1 

(2-njf) : 2 


2 

G ii< f > 


( 2irf ) ” 4 G a (f) 


( 12 ) 


It is important to point out that in equation (12) the dimension of G (f) is supposed 
in. /Hz, thereby requiring that G--(f) have the dimension (in. /sec ) / Hz, the use of 
the tilde (~) serving to distinguish between G^(f) and G-(f), which has the dimension 
g 2 /Hz. In terms of G-(f) 


G u (f) = Y f " 4 G U (f) (13) 

where the numerical value of y is given by 
y = (386. 08858) 2 /(2 tt) 4 , 

2 

it being assumed that the local acceleration due to gravity is 386.08858 (in. /sec ). 

In a manner similar to that of arriving at equation (13), one can argue that 

G^f) = y* f" 2 G u (f) . (14) 

The dimension of G^(f) is (in. /sec) /HZ and y* is given by 

y* = (386.08858/2 -tt) 2 . 

Between G u (f) and G^(f) is the obvious relation 

G u (f) = (2 tt f)‘ 2 G^f) . (15) 


Among other obvious relations are the following: 
< ?Si) - -< £?(Q £-vve 

(U) ^TuT ^TuT * 5 x ’ y,e 


T 5/U (S) = T r//, (s) = T r/,i< s > > F > = x >y-° 


C/U 


£ /u 


G--(f) = |T e/u (2njOr G ii (f) * C = x,y 


(16) 
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(17) 


G qC f ) = (386.08858) 2 |T e/u (2njf) | 2 G~(f) 

G 5 (f) = l T 5/u (27ljf) | 2 G ii (f) ’ 5 = x ’ y>0 

G f (f) = |T Uu (2^f)| 2 G u (f) , 5 = x,y,6 


(18) 

(19) 


The numerical factor was introduced in equation (17) because, as mentioned before, 
the dimension of G~(f) is g 2 /HZ. 

The most efficient sequence of instructions to be executed in computing the mean 
squares and root mean squares of both input and output variables is the following: 

1. Assign a value to f and compute G~(f) in accordance with the expressions 
(defining the curve fit) appearing beneath the hypothetical plot of Figure 3. 

2. Compute and store G^(f) in accordance with equation (14). 

3. Compute and store G u (f) in accordance with equation (15). 

4. Compute and store G”(f) in accordance with equations (16) and (17), 

K = x,y , 9. " 

5. Compute and store G/(f) in accordance with equation (18), £ = x,y,0. 

6. Compute and store G^(f) in accordance with equation (19), £ = x,y,9. 

7. Increase f by Af. 

8. Repeat 1 through 7 until the frequency interval over which G^(f) is pre- 
scribed has been covered. (In this paragraph f ^ and f^ will denote the 
left and right extremes of that interval.) 

9. Via some numerical integration scheme, compute the mean square of £, 

o 

denoted by £ (f^,f^ P er ^ nen l to the interval (f^,f^) in accordance with 



(f l’ f N 



G ? (f) df , 


K = u,u,u,x,x,x,0, e, e,y,y ,y 


( 20 ) 
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10. Extract the square root of £ (f, ,f N ) to get the root mean square (RMS) 
of 5. 


^RMS^W 


= { c 2 (f 1 ,f N )| 1/2 , 5 = u,u,u,x,x,x,y,y,y,'e,6,9 
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In numerically evaluating the integral in equation (20), the author has found 
that the simple trapezoidal rule gives satisfactory results, provided a wise choice of 
Af is made; but, at this writing, can offer no failure proof method for selecting the 
"optimum" value of A f in a given case. Usually, one relies on experience * in deciding 
the value to be assigned to Af. 

To find the mean squares of ii, u, and u it is not necessary to resort to any 
numerical integration scheme since closed expressions are available for their computa- 
tion. From Reference 1 


u2 ( f l» f N ) ~ 1 + b. {f EX,i+l G ii (f EX,i+l* " f EX,i G u (f EX,i )} 


(b.#-l) 


? •‘"(fer) 


, (1 < i < NSEG) (21) 


(b.=-l) 


u (f r f N ) 


“ Y ' 23 b. 1 -! { f EX,i+l G ii^EX,i+P ‘ f EX,i G u (f EX,P } 


(b.#l) 


+ v« 


E , / f EX,i+l\ 

Ci 1^x7”/ ’ 


(1 < i * NSEG) , (22) 


(b i= l) 


U (f l ,f N ) " Y Z b .-3 { f EX,i+l G ii (f EX,i+l ) " f EX,i G u (f EX,i ) } 


(b.#3) 


+ Y 


£ c. In K — ’ 1 — ) , (1 ^ i ^ NSEG) . (23) 

• \ EX,i / 


(b.=3) 


1. A visual examination of the plot of G^(f) could be of some use in deciding whether 
to pronounce a specific value of A f as satisfactory or unsatisfactory. 
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In equations (21), (22), and (23), the symbols . and . + ^ denote, respectively, 
the abscissa of the left extremity and right extremity of the ith straight line segment 
in the log-log plot of G^(f), there being NSEG such segments (see Figure 3), and 
ff = fgx p fjg = fgjj NSEG+1* N °ti ce that equations (22) and (23) have meaning 
only if f^ > f^ >0, and further, that when G~(f) = W(g 2 /HZ) = a constant for 
^ i f < f N , equations (21), (22), and (23) become equations (21)', (22)’, and (23)', 
respectively . 


u 2 (f r f N ) 

u 2 (f r f N ) 


= W(f N -f 1 ) 


(21)' 

= y' W(f" 1 1 - 

f~ *) 

(22)' 

_ yW , f -3 
~ 3 v 1 

f- 3. 
N } 

(23)’ 


While dwelling on "closed expressions," mention should be made of the existence of 
closed expressions for x, 0, 0, 6, y, and y in the very special case wherein G^(f) 
is constant ^ for 0 < f < » . In this case, it is not difficult, with the aid of the table 
of integrals in Reference 2 (see also References 3 and 4, both of which cite Reference 
2) , to show that the mean square of £ , pertinent to the semi-infinite frequency inter- 
val (0,“), is given in closed form by 


? 2 ( 0 ,°°) 


Y*W 


{< B 0 /A 0> (A 2 A 3 - A lV + A 3 <B 1 - 2 B 0 B 2> 

+ A 1 < B 2 - 2 W + < B 3 /A 4> < A 1 A 2 - W } (;) > {' A 0 A 1 
+ A x (A 2 A 3 - A^ 4 ) | , 5 = x,e,e,0,y,y , (24) 


where W denotes the constant value of G^(f) and the numerical factor y* depends 
upon which of the variables £ represents, that is, 

1 1.0 if £ = x 

(386.08858) 2 if £ = 6,0,0,y,y 

2. In the jargon of vibration engineers the base acceleration in this case is termed 
"white noise." 
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The subscript £ on the right brace in the numerator of equation (24) serves to 
indicate that the (K = 0,1,2, 3) are pertinent to the particular £ being dealt with. 

** tr 

The Aj^ (K = 0,1, 2,3,4) are the same for all 5, A K being the coefficient of S in the 
system characteristic polynomial, A(S), defined by the last of equations (9). On 
performing the indicated multiplications in equation (9) and collecting terms one will 
find 


A 0 = K 1 K 2 (L 1 + L 2 )2 

A 1 = (C 1 + C 2 ) (K 1 L 1 2 + K 2 L 2 2) + (K 1 + K 2 } (C 1 D 1 2 + C 2 D 2 2) 
- 2 (C 2 D 2 - CjDj) (K 2 L 2 - K^) 

A 2 = m (K 1 L 1 2 + K 2 L 2 2 ) + I (K x + K 2 ) + (D x + D 2 ) 2 

A 3 = m (CjD^ + C 2 D 2 2 ) + I (C 1 + C 2 ) 

A 4 = ml . 


Pertinent to x, the B^, K = 0,1, 2, 3, are 



B 2 = A 2 - m (K^ 2 + K 2 L 2 2 ) 

B 3 = ’< C 1 + C 2> 

Pertinent to 0, the B^, K = 0,1, 2, 3, are 
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Pertinent to 0, the B^, K = 0,1, 2, 3, are 


B 1 = m ( K 2 L 2 ' K 1 L 1> 
b 2 = 0 

B 3 = 0 

Pertinent to 6, the B K> K = 0,1, 2, 3, are 
B 0 = m (K 2 L 2 - K^) 

B 1 = B 2 = B 3 = 0 * 

Pertinent to y, the B R , K = 0,1, 2, 3, are 

B o = ° 

= m (KjLj 2 + K 2 L 2 2 ) 

B 2 = m CCjDj 2 + C 2 D 2 2 ) 

b 3 = ml . 

Pertinent to y, the B R , K = 0,1, 2, 3, are 


B 1 =m (CjD^ + C 2 D 2 2 ) 

B 2 = mI 
b 3 = ° . 

The structure of the transfer functions relevant to x, x, and y is such as to pre- 
clude use of the referenced list of integrals to find the mean squares of x, x, and y. 

Equations (1) through (23), plus attendant relations (Appendix A), constitute 

the basis for program TRROBM (a mnemonic for "Translational and Rotational Response 

to Base Motion") which has been operational since 1983. A recent revision of the 1983 

version was made so that the program output would include items of importance to the 

author in dealing with a related assignment. Before further comment regarding the 

related assignment is made, the author would like to call attention to Table 1 which 

shows the remarkably close approximations, afforded by equation (20), to the mean 
2 

squares E, (0,°°), E, = x,0,0,6,y,y, whose exact values are determined by equation (24). 
Below the table are the specifications defining the case which was processed by pro- 
gram TRROBM to get the entries in the third column. The coding of program TRROBM 
requires that the input include the items appearing in the left hand column of Table 2. 
Consequently, when certain of the system parameters are "indirectly" specified, as in 
the manner beneath Table 1, one must resort to some preliminary computation in 
accordance with the equations of APPENDIX B to determine the numerical values of 

^ 2 ’ ^p ^2* B p ^2* f * 

Not shown in the list of input items in Table 2 are other input items which are 
"implied" by the presence of G^(f) in that list and by the expressions for the curve 
fit parameters under the diagram of Figure 3. These items include NSEG, f.™ . 

(i = 1, . . . ,NSEG+1), NCORN , GCORN , and ADB(i,i+l) . i = 1,...,NSEG, all essential in 
the computation of G-(f) for a given value of f. In program TRROBM the f^ x ^ and 
ADB(i, i+1) are embedded in the one- dimensional arrays identified by the FORTRAN 
symbols FEX and DELDB , respectively. By mere inspection of the PSD specification, 
one has immediately the input data designated NSEG, FEX, DELDB, NCORN, and 
GCORN. Pertinent to the sample PSD of Figure 2 these items are 
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NSEG = 7 


f EX,l 


20. 


’ +6. 

f EX,2 


30. 


0 

f EX,3 


120. 


+6. 

f EX,4 

= 

210. 

(HZ), DELDB = 

0 

f EX,5 


400. 


+9. 

f EX,6 


480. 


0 . 

f EX,7 


900. 


_ -!2. 

1 

00 

>< 

< 4-1 

1 


2000. 




NCORN = 1 , GCORN = 0.15 (g 2 /HZ) 


(DB /OCTAVE) 


The choice of the combination NCORN = 1, GCORN = 0.15 was but one of several 
available. The admissible combinations of NCORN and GCORN in this case are shown 
in the following table. 


NCORN 

1 

2 

3 

4 

5 

6 

7 

8 

GCORN 

0.15 

0.32 

0.32 

1.0 

1.0 

1.7 

1.7 

0.075 


It is evident from Table 2 that the entries in the second column of Table 1 are 
not to be found among the items output by program TRROBM. Instead, they are the 
output of a smaller program, an auxiliary to TRROBM (aptly named program AUXRBM), 
which was coded only recently, in July 1986. It was long after the author had 
developed program TRROBM and two similar programs that he learned, through 
browsing the literature (References 2, 3, and 4 in particular), of the existence of 
the table of integrals which served as a guide in writing equation (24) upon which 
program AUXRBM is based. Neither was he aware, until he surveyed the literature, 
that much of the work he had done in developing the two programs described in the 
footnote had already been done years ago. 

3. Program RESPBM (response to base motion) treats the single d.o.f. mass-spring- 
damper system excited by the random vibratory motion of the base whose accelera- 
tion PSD is prescribed as in this paper. Program RESBM2 deals with the randomly 
base driven 2 mass- 2 d.o.f. system, the two d.o.f.'s being translational. 
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Program AUXRBM, a sample output of which is given in Appendix C, provided 

the numerical data necessary to the construction of the families of curves in Figures 

4 through 11. The data could have been generated by program TRROBM but at a 

greater cost of computer time, not to mention the slight inaccuracies in the data due 

to the necessity of restricting the mean square computation to a finite frequency 

interval whose left extremity must be positive. The use of the word "inaccuracies" 

tends to unjustly discredit program TRROBM. In defense of TRROBM the author 

should point out that even in those cases wherein G^(f) is constant, which is the only 

kind of case to which AUXRBM is applicable, there is hardly a discernible difference 

between the plots 4 5 6 of RMS's made from the output of TRROBM and those made from 

the output of AUXRBM (after the numerical values have been rounded to at most 

three significant digits and plotting is done using the same scales for both sets of 

output). The author has made this assertion on the assumption that, in processing 

a case by TRROBM, a wise choice of Af is made, and further, that a sufficiently 
£ 

wide frequency interval is used in the mean square computation. As support of his 
assertion, the author invites the reader to compare the values of 0 R |yjg» e R] yjg 0 rmS’ 
y RMg , -^RMS’ an< * *RMS ^ oun(i amon g the items of the sample TRROBM output in 
Appendix D with the appropriate encircled values or inset tabular entries of Figures 
6 through 11. 

The source of the RMS's in Figures 12, 13, and 14 was program TRROBM. In 

each of these figures the prevailing conditions are the same as those pertinent to the 

encircled points of Figure 6. The previously cited tables of integrals, and hence, 

"2 1 TTo 

program AUXRBM, were of no utility in the computation of x , x , and y because, 
as mentioned in a previous paragraph, the rational functions 


x/ii 


(s > = 3 T */u (s) 


T 


x/ii 


(s) 



T x/u (s) 


9 


4. For the purpose of comparing the behavior of the plotted function, as depicted 
by the solid curves , with its behavior under slightly different conditions , some 
of the figures have either an inset table of values or encircled values of the 
function corresponding to the changes in system parameters. 

5. Plots of RMS's versus r^ (holding r L constant). 

6. Confining the mean square computation to the interval (1., 2f 2c ) results in 

excellent approximations. ^ 
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and 


T 


y/ii (s > = T y/u 


(s) 


do not have the requisite structure. The first thing one will notice about these figures 
is that no attempt has been made to draw a "best fit" curve through any of the 
several sets of points, the reason being an insufficient number of points to accurately 
determine the behavior of the variable plotted. 

All of the programs described in this paper are coded in FORTRAN V for a 
punched card machine (one of the UNIVAC series in particular). However, one with 
the expertise can translate the FORTRAN language into that of another computer. 

A fellow employee here at MSFC has, in fact, already effected the translation of pro- 
gram RESPBM, described in one of the footnotes, into TEKTRONIX language (models 
4051, 4052, and 4054). 

When the author was approached by his supervisor with questions about the 

g 

mean values of system kinetic energy, potential energy and energy dissipated, his 
first thought was of the system coordinate velocities whose mean squares are not a 

a 

part of the output of the 1983 version of TRR OBM . It was the need for the mean 

o 

squares of the coordinate velocities, as well as yp (whose need will become apparent 
later), that prompted the 1986 revision of the 1983 version of TRROBM. In response 
to the questions asked, the author has developed the following expressions for the 
mean values of the kinetic energy (KE) and potential energy (PE), the symbol E in 
E(£) being the familiar expectation operator or mean value operator. 


E(K.E.) - \ (m Xp Mg + I 9 r MS ) 


E(P.E . ) = 


m oj 


hx y 2 


I w 


RMS 


2 0 RMS 2 L CK 2 L 2 
P 


W <yp, 


RMS 


2 

y RMS 


' L P 9 RMS ) + 2 (K 1 6 ST,1 + K 2 6 ST, 2* 


7. Pat Lewallen, ED 24. 

8. See Figure 1. 

9. When the work which culminated in the 1983 version of TRROBM was done, interest 
was primarily in accelerations and displacements. 
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y p, and Lp are given elsewhere but are repeated here. 


2 2 

The definitions of u nx> a> nQ , 

“nx = < K 1 + V / m 

“na = < K 1 L 1 2 + K 1 L 2 2) ' 1 


Yp = displacement of arbitrary point P (not the CM) relative to the base 

Lp = lateral distance of point P from the CM, positive or negative according 
as point P is right or left of the cm 

The expression for E(PE) was derived on the assumption that the mean value of u, 
the base displacement, is zero, and that the zero level for gravitational potential 
energy is the static equilibrium level of the CM. Assuming further that = Dg = D 

and C. = C 9 = C, it is not difficult to show that the mean value of the rate at which 

1 * • 2 2 '2 
energy is dissipated through viscous damping is 2C(y RM g + D e RMS ^ * Develo P ment 

of the expressions for the mean values in this paragraph was the "related assignment" 

alluded to earlier. 

~o 

Attention is now called to Figure 4 wherein the symbol {x ( 0 , «> ) ^ q = q denotes 
the mean square of x (for 0 < f < °°) when component rotation has been suppressed 
entirely by enforcing the relations KjL^ - KgL 2 = 0 and CjD^ - CgDg = 0 so that 6 
is identically zero (provided 0 and 0 are initially zero). A cursory examination of 
this family of curves reveals that suppressing rotation merely serves to increase the 
mean square of x (otherwise, the plotted mean square ratios would be greater than 
one) . 


It is evident in Figure 5 that for some combinations of r^. and r^, which admit 
rotation, the mean square of y is larger than it is when there is no rotation while 
for other combinations it is smaller. 

Figures 6 through 11 show whether the imposition of the conditions (D^ = D 2 , 
4 L. , i = 1,2} instead of {D. = L. , i = 1,2}, other conditions being the same , 
results in an increase or decrease in the RMS of the response variable in question. 
These figures, when complemented by the output of TRROBM plotted in Figures 12, 
13, and 14, provide the RMS's of the system coordinates and their first two time 
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derivatives for both the (x,0) and (y,0) descriptions of system configuration. This 
collection of figures does not represent an exhaustive parameter study, but is 
exemplary of parameter studies made possible by program TRROBM (with or without 
the support of its auxiliary AUXRBM, which is of limited application). 

At least one paragraph should be devoted to stability, if only to go so far as 

to write the conditions (on the system parameters) whose satisfaction guarantees 

system stability. Such conditions are indirectly realized by conditions on the coeffi- 

4 

cients of the system characteristic polynomial, A(S) = £ A.S 1 [see equations (9)] , 

i=0 1 

those conditions being available via the Routh-Hurwitz criterion. The Routh-Hurwitz 
array pertinent to A(S) is the following: 


ROW 

A 4 

A 2 

A o 

0 

1 

A 3 

A 1 


2 

A 3 A 2 " A 4 A 1 
A 3 

A o 


3 

A i( A 2A 3 - Al A 4 ) - A 0 A 3 2 
A 3 A 2 " A 4 A 1 



4 

A 0 




By the expressions defining them, A^, A^, and Aq are intrinsically positive. Hence, 
by the Routh-Hurwitz criterion for stability, system stability is assured if the other 
two elements in the first column of the array are positive, or, equivalently, if both 
of the inequalities 

A 3 A 2 > A 4 A 1 ’ A l^ A 2 A 3 ' A 1 A 4^ > A 0 A 3 2 ’ 

are satisfied. The problem of assessing the "degree" of stability will not be addressed 
in this paper. The technique for handling the situation wherein a left column element 
is zero will be found in the literature. 


15 




The author has given some thought to models other than that of Figure 1, 
those of Figure 15 in particular, which, in certain cases, could be more "credible" or 
"plausible" models. Though there is slight difference between the "appearances" of 
the models in Figures 1 and 15-a, that difference due, obviously to the elastically 
supported dampers in Figure 15-a, there is a marked difference between the respective 
mathematical descriptions of model motion. While the differential equations governing 
the motion of the model in Figure 1 are of second order, those determining the motion 
of that in Figure 15-a are of third order. 

Pertinent to Figure 15-a, the author has derived the following equations. 


mCy+ii) = - (K x + K 2 ) y + (K^ - K 2 L 2 ) 6 - (y - D^) - C 2 (y + DgO) 


1 1 8 - m D 2 (y + u) + [K 2 L 2 - - D 2 (K x + K 2 )] y 


+ [KjLjl 2 + k 2 l 2 2 + D 2 (K 1 L 1 - k 2 l 2 )] ej 


- 1*0- m D 1 (y + u) - [K 2 L 2 - K 1 L 1 + D 1 (Kj. + Kg)] y 


- [KjLj 2 + K 2 L 2 2 - D x (K 1 L 1 - K 2 L 2 )] 6 } (25) 


I 0 = (K 1 L 1 - K 2 L 2 ) y - (K 1 L 1 2 + K 2 L 2 2 ) 6 + (y - D^) - C 2 D 2 (y + Dg9) 


C 1 D 1 

' — lie - m D 2 (y + u) + [K 2 L 2 - K 1 L 1 - D 2 (Kj. + Kg)] y 

Ki(Di + D 2 ) ^ 

+ + K 2 L 2 2 + D 2 (KjLj - K 2 L 2 )] 0} 

C,D ( ... 

+ j-ie-mD 1 (y + u) - [K 2 L 2 - K 1 L 1 + Dj (Kj + Kg)] y 

k 2 ( d i + d 2 ) 

- [KjLj 2 + K 2 L 2 2 - Dj (KjLj - K 2 L 2 )I 6 1 . (26) 


k i<Di + d 2 ) 


k 2 (Di+d 2 ) 
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Notice that if one makes the substitution x = y + u and allows and K 2 to become 
infinite, equations (25) and (26) become, in view of relations (3), the equations of 
motion of the model in Figure 1 pertinent to the (x,0) description of system configura- 
tion, the equivalent of the matrix equation (4); or, if one merely permits K.^ and K 2 
to approach infinity, equations (25) and (26) become the equivalent of equation (5). 

Considerable simplification of equations (25) and (26) is realized in the special 
case wherein = K 2 = K, = Cg = C, = K 2 = K. In that case they read as 
equations (27) and (28). 


H) i 


— y + m y + 2C [ 1 + ^ |y + 2Ky + C 
K V K 


-mC 


D 2 - D 1 + § < L 2 - L l>] 9 


+ K (L 2 - L^) 0 = ~ u - m u 

K 


(27) 


IC 

K 


e + I 0 + c 


+ c 


[ d i 2 + D 

h' 


2 * /t 2 + y 2. 

2 + K (Ll L 2 ^ 


j 0 + K (L x 


2 + L 2 2) 8 


K 


>] y 


D i + ^ (L 2 ■ L i ) I y + K (L 2 - L i> y = 0 • 

K 


(28) 


Having written the equations of motion, the next step toward a mean square computa- 
tion is the deduction of the relevant transfer functions. Pertinent to the system 
comprised of equations (27) and (28), it is easily deduced that 


5 y,U< S > = L h S ' ' L l i Sl 

i=0 i=0 


5 e/u (s > = Z ?i s ‘ ' E X i s ‘ 
i=0 i=0 


(29) 


where 
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A 0 = K 2 <Lj + L 2 ) 2 


A 1 = 2 K C 


[ D i 2 * D 


2 2 - (D 2 - Dj) CL 2 - L t ) + Lj 2 + L 2 2 + 


| < L 1 + L 2 )2 ] 


A 2 = m K (Lj 2 + L 2 2 ) + 2 K I + 2 C 2 


- C' 


H)[ 

[ D 2 - D 1 + | (L 2 - L l>] 


D l 2 + D 2 2 + | (L 1 2 + L 2 2 > 
K 


(30) 


A mKC , T 2 ^ T 2^ j. 2 K C I m ^ 
A 0 = — (L-^ + L 2 ) + — + m C 


K 


K 


[ D i 2 + D 


2 2 + | <L 1 2 + L 2 2 > 




A 4 = 


m C' 
K 


[ D l 2 + °2 2 + | (L 1 2 + L 2 2) j + m 1 + ( X + |) 


A 5 = 2 m I C / K 


A g = m I C 2 / K 2 


3 0 = &i = 0 


?2 = - m K (L x 2 + L 2 2 ) 


^ "2 m C K f T 2 ^ T 2 ^ p ^ 2 x ^ 2 ^ 

H = (L x + L 2 ) - m C (D x + D 2 ) 


g 4 = 


-m C 
K 


K 

2 


[ D i 2 + D 


2 , K 2 , r 2% | t 

2 + ^ (L 1 + L 2 ) I " m 1 

K 


i- 


( 31 ) 


e 5 = - 2 m C I / K 

^ 9 ^9 

Ug = -m I C / K Z 
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Y 0 = Yi = 0 


Y 2 = m K (L 2 - Lj^) 

Y 3 = — ^ - C ~ (L 2 - L x ) + m C (D 2 - D x ) 

K 


y 4 = 


m n 2 
m C 

K 



D 1 + | (L 2 - L l>' 


(32) 


Imposition of the additional conditions and D 2 = as in Figure 15-b, results 

in a simplification of equations (27) through (32). Further simplification is possible 
by setting K = K. 

Upon examining the structure of the transfer functions in equation (29), along 
with that of the five Laplace transform ratios 


= ±1 T , (S) 

^f(ii) S 2 y/u 


^(y) _ i 5; 

5J(ii) S y/u 


(S) 


2(6) . 

dC(vi) " 


^ 5 e/u (s > 


^(9) 

St (ii) 


E T e/u (S) 


jj (j) _ ™ 

5Tu) " T e/u tS) 


it_is evident, in light of the table of integrals previously cited, that the mean squares 
— 2 • • 

£ , K = y, y» 9, 9, 0, are expressible in closed form when G~(f) is constant for 

U 2 

0 i f < *, In the notation of this paper, the closed expression for £ is 
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, ( 33 ) 


« 2 = TTZ ®5 2 n 0 + n l (S 4 2 - 2 W + n 2 <S 3 2 - 2 B 2 B 4 + 2 W 
6 

+ n 3 <V - 2 S 1 S 3 + 2 5 0 5 4> + n 4 «1 2 - 2 W + n 5 S o 2 > (?) 

K = y,y,e,e,e 

the significance of the factors y* and W being the same as in equation (24). The 
constant W, incidentally, is herein assumed to have the dimension g /HZ. The sub- 
script £ on the right brace in equation ( 33) , as one should expect , indicates that 
the B., i = 0, 1, ..., 5, are associated with the £ in question. The A., i = 0, 1, 

..., 6, are given by equation (30) while Ag and the n., i = 0, 1 5, are defined 

by 

n 0 = Y~ ( ^4 n l ‘ A 2 n 2 + A 0 n 3> 

o 

n l = Aq (Ag - A^Ag) + (AjA 4 - AgAg) 
n 2 - A0A3A5 + A^ (A^Ag - AgAg) 
n 3 = A5 (AqAj - A X A 4 ) + Aj^AgAg 

n 4 = — (A 2 n 3 - A 4 n 2 + A g n^ 

A 0 

n 5 = (Ag n 4 - A 4 n 3 + A g n 2 ) 

A 0 

Ag = Aq (Aj^ n g - Ag n 4 + Ag ng) 

Pertinent to y the B . , i = 0, 1, ..., 5, are given by 

B . = 0. „ , i = 0,1, 2, 3, 4 [the 0., i = 0,1,. ..,6, are given by equations (31)] 



Pertinent to y. 


B i = ^i+l ’ i = 0,1,..., 5 
• » 

Pertinent to 0, 

B- = y. , i = 0,1, 2, 3, 4 [the y^, i = 0,1,..., 4, are given by equations (32)] 

i 5 = 0 . 

Pertinent to 0, 

Bj = » i = 0,1, 2, 3 

B 4 = B 5 = 0 ‘ 

Pertinent to 0, 

B i = Y i+2 * 1 = 0,1,2 

S 3 = 5 4 = 5 5 = 0 • 

To date, no attempt has been made to code a program based on equations (25) 
and (26) or any of their simplified forms. Programming, on the part of the author, 
has been pursued only so far as programs TRROBM and AUXRBM, which mark the 
culmination of the author's effort in this area. 
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TABLE 1 


£ £ 2 (0,°°) £ 2 (1., 2f 2 C ) Percent Error 

x 10 3 (0.78567859) (g 2 ) 10 3 (0.78498616) (g 2 ) 0.088 

0 10 4 (0.20118432) (rad. /sec 2 ) 2 10 4 (0.20111040) (rad. /sec 2 ) 2 0.037 

0 10' 1 (0.50960584) (rad. /sec) 2 10 _1 (0.51124543) (rad. /sec) 2 0.32 

0 10' 4 (0.12217110) (rad.) 2 10 _4 (0.12244709) (rad.) 2 0.23 

y 10 3 (0.29642678) (in. /sec) 2 10 3 (0.29412964) (in. /sec) 2 0.77 

y 10' 3 (0.75068497) (in.) 2 10" 3 (0.75001788) (in.) 2 0.089 

VIA EQ. (24) VIA EQ. (20) 


SPECIFICATIONS: 

G u (f) = 0.1 (g 2 /HZ) , 1 < f < 200.0 

K 1 /K 2 = 1., Cj/Cg = 1., Dj/Dg =1., ?x = c 0 = 0.01, f nx = 100. (HZ), 


r f = 0.1, r L = 2/3 


m = 1.0 


, lb.*sec 2 
v in. 


p = 5. (in.) 


VALUES OF I, K 1 , K„, C., C 9 , D. , D 9 , L. , L 9 (ENFORCED BY THE 
SPECIFICATIONS): 

I = 25. (lb. *sec 2 *in. ) 

K 1 = K 2 = 10 6 (0.19739209) (lb/in.) 

C 1 = C 2 = 10 (0.62831853) (lb/(in. /sec)) 

= D 9 = 10 (0. 15811388) (in.) 

L t = 0.39223227 


L 2 = 0.58834841 (in.) 




TABLE 2 


PROGRAM TRROBM 


Program Input 

Program Output* 

G a (f) 

G,( f ) 


Af 


► £ = u,u,u,x,x,x,’e\e,e,y,y,y,y p ,y p ,y p 


? rms (f l ,f N ) 



J 


C 1 

f ic * i = 1 ’ 2 

C 2 

Modal Column Corresponding to fj * i = 1,2 

D 1 

f nx 


D 2 

f ne 


K 1 

l T «/u (2 "i f) l 2 

, E, = x, 6,y ,y p 

K 2 


L 1 



L 2 



L P 



m 



I 




* A print of all tabulated functions of frequency is optional 
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plot er sample input psi> (inset) 



FIGURE 2 
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/ _ J _ ADB&i+i) 

P* ™ r /* ^ — * 
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/OA? 2 
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DEFINITIONS 


c. 

1 

Damping coefficient [lb /(in. /see)] of viscous damper i, i = 1, 

D. 

1 

Lateral distance (in.) between CM and point at which damper 
attached to model, i = 1,2 (D^ > 0) 

r D 

D l /D 2 

f 

Frequency (HZ), fj i f < f N 

f. 

1C 

Frequency (HZ) of undamped coupled natural mode i, i = 1,2 

f nx 

Frequency (HZ) of undamped uncoupled translational mode 

f ne 

Frequency (HZ) of undamped uncoupled rotational mode 

g 

o 

Acceleration due to gravity (in. /sec ) 

G ii ( f ) 

PSD or ii (g 2 /HZ) 

G-(f) 

PSD of u [(in. /sec) 2 /HZ] 

G u (0 

PSD of u (in. 2 /HZ) 

G x< f > 

PSD of x (g 2 /HZ) 

c x (f) 

PSD of x [(in. /sec) 2 /HZ] 

G x (f) 

PSD of x (in. 2 /HZ) 

G.-(f) 

PSD of 0 [(rad. /sec 2 ) 2 /HZ] 

G i(f) 

PSD of 0 [(rad. /sec) 2 /HZ] 

G 0 (f) 

PSD of 6 (rad. 2 /HZ) 

G y(f) 

PSD of y (g 2 /HZ) 

G y (f > 

PSD of y [(in. /sec) 2 /HZ] 

G y tf) 

PSD of y (in. 2 /HZ) 

G- (0 

*P 

PSD of y p (g 2 /HZ) 
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6 <f r f N > 

e rms (f l ,f N ) 


9 (f r %) 


0 rms^l’*N^ 


n f r f N ) 


0 rms (f l’ f N> 


y 

y 

y 


y (f r f N ) 


y rms^l ’^N ^ 


y ( f x» f N ) 


y rms (f l ,f N ) 


y < f i» f N ) 


y rms (f l> f N> 


Angular displacement (rad.) of model from static equilibrium 
orientation 

Angular velocity (rad. /sec) of model 

2 

Angular acceleration (rad. /sec ) of model 

2 

Mean square value (rad. ) of 6 in the interval (f^>f^) 

Root mean square (rad.) of 0 on the interval (fpfjj) 

Mean square value, (rad. /sec) , of 0 on the interval 

Root mean square (rad. /sec) of 0 on the interval (fpfjj) 

2 2 

Mean square value [(rad. /sec ) ] of 0 on the interval 

2 

Root mean square (rad. /sec ) of 0 on the interval (fpf N ) 

Displacement (in.) of model CM relative to the base 

Velocity (in. /sec) of model CM relative to the base 

2 

Acceleration (in. /sec ) of model CM relative to the base 

2 

Mean square value (in. ) of y on the interval ( f ^ , f ^ ) 

Root mean square (in.) of y on the interval 

2 • 

Mean square value, (in. /sec) , of y on the interval (fpf N ) 

Root mean square (in. /sec) of y on the interval (fpf N ) 

2 

Mean square value (g ) of y on interval (f^f^) 

Root mean square (g) of y on the interval 
Displacement (in.) of point P relative to the base 
Velocity (in. /sec) of point P relative to the base 


y P 

y p 2(f l ,f N ) 

y P,rms (f l ,f N ) 

y P (f l ,f N ) 

y P,rms (f l’ f N ) 

y p 2(f l ,f N ) 

y P,rms (f l ,f N ) 

6 ST,i 

0 ). 

1C 

Af 


U 

u 

ii 


vAf r f N ) 


U rms (f l ,f N ) 


« (f r f N ) 


u rms (f l’ f N ) 


u"(f r f N ) 


u rms (f l ,f N ) 


2 

Acceleration (in. /sec ) of point P relative to the base 

2 

Mean square value (in. ) of yp on the interval 
Root mean square (in.) of yp on the interval 

2 

Mean square value, (in. /sec) , of y p on the interval ( f x , f N ) 

Root mean square (in. /sec) of y p on the interval (fpf N ) 

2 

Mean square value (g ) of y p on the interval (f^.f^j) 

Root mean square (g) of y p on the interval 

Static deflection (in.) of spring i, i = 1,2 (considered positive 
despite my sign convention) 

2irf ic (rad. /sec), i = 1,2 

Both the print step (HZ) and the frequency increment (HZ) used 
in the numerical evaluation of the definite integrals defining the 
mean squares 

Displacement (in.) of the base from its static equilibrium position 
Base velocity (in. /sec) 

2 

Base acceleration (in. /sec ) 

2 

Mean square value (in. ) of u on the interval (*i»f N ) 

Root mean square value (in.) of u on the interval (fpf N ) 

Mean square value (in./sec)^ of u on the interval 
Root mean square value (in. /sec) of u on the interval 

2 

Mean square value (g ) of u on the interval (fpf N ) 

Root mean square value (g) of u on the interval 

Displacement (in.) of model CM from its static equilibrium position 
(x is an absolute displacement) 

Absolute velocity (in. /sec) of model CM 
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X (fpf N ) 


X rms (f l ,f N ) 


X ^1’%^ 


X (f,,^.) 
rms v 1 N' 


x ‘(f r f N ) 


x rms (f l ,f N ) 


|T ?/U (2rfi)| 2 


T e/U (2rfj) 


K. 

l 

K. 

l 

L. 

l 


r 


L 


m 


P 


L 


P 


2 

Absolute acceleration (in. /sec ) of model CM 

2 

Mean square value (in. ) of x on the interval (fpf^) 

Root mean square (in.) of x on the interval (fpf N ) 

2 

Mean square value (in. /sec) of x on the interval (f^f^) 

Root mean square value (in. /sec) of x on the interval 

Mean square value (g^) of x on the interval 

Root mean square (g) of x on the interval (f^,f^) 

The complex variable of the laplace transformation (rad. /sec) 

Square of the magnitude of the frequency response function between 
£ and u (dimensionless), £ = x.y.yp (j = AT! ) 

Square of the magnitude of the frequency response function between 
0 and u (rad. /in. )2 

Moment of inertia about an axis through the CM and perpendicular 
to the plane of motion (lb*sec^*in. ) 

Stiffness (lb /in.) of linear spring i, i = 1,2 

Stiffness (lb /in.) of elastic support of damper i in alternate model 
(Fig. 15), i = 1,2 

Lateral distance (in.) between model CM and the point at which 
spring i is attached to the model, i = 1,2 

L l /L 2 

2 

Mass of component [lb/(in./sec )] 

Radius of gyration (in.) 

Lateral distance (in.) between model CM and point P (L p > 0 or 
Lp < 0 according as point P is right or left of the CM) r 

f /f 
n6 nx 

Fraction of critical damping associated with translation 
Fraction of critical damping associated with rotation 
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NSEG 

f EX,i 

NCORN 

GCORN 

ADB(i,i+l) 


Parameters appearing in the analytical representation of G”(f) on the 
interval (f EX>i > f EX,i+P’ i = !»**** NSEG 

Number of straight line segments in the log-log plot of the prescribed 
base acceleration PSD 

Abscissa (HZ) of the i'th "corner" point in the log-log plot of G-j(f), 
i = 1, . . .NSEG+1 

A positive integer specifying that "corner" point of the straight line 
segment representation of the input base acceleration PSD on log-log 
graph paper at which the value of the input base acceleration PSD 
is given 

o 

The value (g /HZ) of the input base acceleration PSD at the "corner" 
point specified by the integer NCORN 

Rate of change, in decibels /octave , of the input base acceleration 
PSD as the frequency f varies from f £X . to f EX . +1> i = 1....NSEG 
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APPENDIX A 


SUBSIDIARY RELATIONS 


R 2 . a 2 

|T x/u (2 1 ,jO| 2 = X 


{Re(D)} 2 + U#m(D)} 2 


R x = K 1 K 2 (L 1 + L 2 )2 " [I(K l + K 2 } + c i C 2 < D i + d 2 )2] (2irf)2 


*4c = " 1 (C 1 + C 2 ) (27rf) 3 + (2rf) t(C x + C 2 ) (K x L x 2 + K 2 L/) 

+ (K x + K 2 ) (C 1 D 1 2 + C 2 D 2 2 ) + 2 (K 1 L 1 - K 2 L 2 ) (C 2 D 2 - C 1 D 1 )] 


R 0 2 + ^e 2 


l T e/u (2 ^r = 2 a 2 

/U {Re(D) i + {\f/ m(D) } 2 


R 0 = - m (K 2 L 2 - K 1 L 1 ) ( 2 uf ) ' 


= - m (C 2 D 2 - C 1 D 1 ) (2wf) : 


Re(D) = m I (2Trf) 4 - [m (K^ 2 + KgLg 2 ) + I (K x + Kg) + (D t + D 


+ K lK 2 «*i + V 


^m(D) = f ( C x + C 2 ) (K^j 2 + K 2 L 2 2 ) + (K 1 + Kg) (C ^ 2 + CgD^) 


+ 2 (K 1 L 1 - K 2 L 2 ) (C 2 D 2 - CjDj)] ( 2itf ) - [m (C^ 2 + CgDg 2 


+ 1 (C i + C 2 )] (2rf )' 


2 ) 2 ] (2TTf) 
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FOUR ADMISSIBLE EXPRESSIONS FOR THE 
UNDAMPED MODAL MATRIX (NON -NORMALIZED) 
(K L + K 2 L 9 PRESUMED) 
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APPENDIX B 


ON THE CONSEQUENCES OF CERTAIN SPECIFICATIONS 

When the coupling coefficients vanish simultaneously, that is, when - 

KgLg = ® an< * ^1^1 " ^2^2 = e Q ua ti° ns of motion, (1) and (2), assume their 

uncoupled forms (B-l) and (B-2). 

mx + (C : + C 2 ) x + (K x + K 2 ) x = ( C 1 + C 2 ) u + (K x + Kg) u (B-l) 

16 + (C 1 D 1 2 + C 2 D 2 2 ) 9 + (KjL^ 2 + K 2 L 2 2 )0 = 0 (B-2) 


Upon inspecting equations (B-l) and (B-2), it is easily seen that the uncoupled 
undamped natural frequencies and damping ratios satisfy 


0) 


2 _ K 1 + K 2 


nx 


m 


2 ? x “nx 


C 1 + C 2 
m 


“ne = r <K 1 L 1 2 + K 2 L 2 2 ) > 2 h “ne = T (C 1 D X 2 + C 2 D 2 2 > • 


If one imposes the conditions 


K 1 K 2 ’ C 1 C 2 ’ D 1 D 2 ’ Sc c 0 » 


and assigns values to m, p, w , x, , r~, r T , where 

nx x i Li 


r f = "ne'V * W f nx and r L * L l /L 2 


then some simple algebraic manipulation (bearing in mind the definition p = A/m) 
will show that the numerical values of D., C^, Lj, and K., i = 1,2, are determined by 


D. = p t/f^T , i = 1,2 


C i = m Sc “nx ’ 1 = X ’ 2 
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If the condition D.^ = D 2 is replaced by Dj = L., i = 1,2, the numerical values 

of G--(f), m, p, a) , l , r», and r T being prescribed as before, the expressions for 
u. nx x i Li 

I, Kp K 2> C 1 , C 2> Lp and L 2 are the same as before, but the equality of and 
L., i = 1,2, requires that ^ In this case one will find 


T x/u (s > = Z B K sK/ Z a k S ' 


, “nx r f 2 (1 + r L> 2 
0 = 2(1+ r L 2 ) 


. 2 s x “nx r f 2 (1 + r l/ 

A 1 = 2 “ “ 

1 + r L 


2 r 2 r f 2 (1 + r ) 2 
1 + r 2 + X ? * ±— 

1 + r L 


» A, = 2 c u> (1 + r f ) , 


x nx' 


A 0 ’ 


B 1 A 1 * B 2 " A 2 " “nx r f 2 ’ B 3 2 c x “nx 


T e/u <s) = Z % s K / Z a k s 1 


^ 0 ~ ' 1 “ 0 ’ ' 2 


IV w 
t nx 


1 - r r 



y 3 


2 r f c a) 
f s x nx 


1 - r r 


/ 2 (1 + r L 2 ) 


4 4 

T y/u^ 53 e K S ^ 53 A k S 

K=0 K=0 


K 


*0 = Pi = 0 ’ *2 = 


r f 2 w nx ’ e 3 " " 2 r f 2 c x w nx ’ 3 4 _1 


If neither the condition D 1 = D 0 nor the condition D. = L. , i = 1,2, is imposed, 

LA 1 1 

while enforcing the conditions = K 2> = C^> and prescribing the numerical 

values of G--(f), m, p, a> , t, , c, a , r-, r T , and r~ = D./D,, there is still no change 
u nx X U I Lt D LA 

in the expressions for I, K^, K 2> C^, and C 2> but and D 2 will be determined by 


°2 = p 


2 r f 


Sc + V > 


1/2 


D 1 r D D 2 


In this special case, the transfer functions relevant to x, 0, and y are 

T x/u< s > ■ E V ' i V ®“ 

K=0 K=0 

A k ' = A^. of the preceding paragraph for K = 0,3,4 


A, = 2 c u> r., 
1 nx f 


+ r - (1 r D> (1 - r L> 


L C X 


h r r i 

(1 + r D 2 ) (1 ♦ r L 2 )J 


1/2 


A„ = o> 

2 nx 


1 + r f “ + 


2 . 2 r f Sc «e (1 + r D>‘ 


1 + r 


D 


B k ' = of the preceding paragraph for K = 0,3 
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It is worthy of note, insofar as economy of computer time is concerned, that the 
moduli of the transfer functions of this and the preceding paragraph are unchanged 
(for a specific value of the complex variable s) if both r L and r^ are replaced by 
their reciprocals. The same can be said of the transfer functions defined by equa- 
tions (6), (7), and (8) when their numerator and denominator coefficients are such 
as to meet specifications similar to those beneath Table 1. 
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APPENDIX C 


SAMPLE OUTPUT OF PROGRAM AUXRBM 

On the following- pages of this Appendix are two sets of mean squai’es and RMS's 
of x, 0, 0, 0, y, and y as found by program AUXRBM in two computer runs, there 
being seven cases processed in each run. In both runs r L assumes the values 2/3, 
1.0, 1.05, 1.2, 1.5, 2.0, and 10, in turn, and 

r f = °' 1 ’ f nx = 100 ‘ (HZ) ’ Sc = °* 01 ’ 

Ki = K 2 , Ci = C 2 , p = 5. (in.) , m = 1.0 (lb*sec 2 /in. ) , 

G”(f) =0.1 (g 2 /HZ) for 0 < f < » . 

In one run the conditions = ^2 and (i = 1,2) apply, while in the other 

D. = L. (i = 1,2). 

Observe that in each run the output mean squares in case 5 duplicate those of 
case 1, that being due to the fact that the values of r L and r D in case 5 are the 
reciprocals of those in case 1. See Appendix B (last two sentences). 

Notice also that the numerical results in case 2 of each run verify that 0 is 
identically zero (assuming zero initial conditions) when the relations C^D.^ = C 2 D 2 and 

K 1 L 1 = K 2 L 2 hold - 

The following table will serve to define the FORTRAN symbols appearing on the 
AUXRBM printout. 
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APPENDIX D 


SAMPLE OUTPUT OF PROGRAM TRROBM 


The specifications G~(f) = 0.1 (g 2 / HZ), m = 1.0 (lb*sec 2 /in. ) , p = 5. (in.), 
f nx = 100. (HZ), C x = C Q = 0.01, K 1 = K 2 , C x = C 2 , D x = D 2 , D. 4 L. (i = 1,2), 
r^ = 2.0, 1. < f < 400., Af = 0.5, and 


r 




1.05, (Case 1) 
1.2 , (Case 2) 
1.5 , (Case 3) 
2.0 , (Case 4) 
10. , (Case 5) 


led to the numerical values of I, K^, K 2> L^, L 2> C^, C 2 , D^, D 2> and other items 
(with the exception of L ) essential to the mean square computation, shown on the 

IT 

input print which precedes the output print of program TRROBM. In each of the 
five cases processed by TRROBM, the frequency interval (1., 400.) HZ was slightly 
less in length than the recommended interval (1., 2f 2c ) HZ, but the approximations 
to the mean squares and RMS's were surprisingly good. The reader should compare 

0 RMS’ 0 RMS* 0 RMS ’ y RMS’ ^RMS’ and *RMS tound in ,he ° utput print * with the 
encircled values and inset tabular values, corresponding to r^ = 2.0, in Figures 6 

through 11. 

The reader is due an explanation of the items appearing on the last page of the 
output print for each case. As it pertains to matrices, the word adjoint has its 
usual meaning, that is, the adjoint of a matrix is the transpose of the associated 
matrix of cofactors. The 2x2 matrix identified as "ADJOINT CORRESPONDING TO 
OMEGA1C" is merely the adjoint of the characteristic matrix, to be defined subse- 
quently, when the elements of the characteristic matrix are evaluated at to = to^£. 

A similar statement applies to the matrix identified as "ADJOINT CORRESPONDING TO 
OMEGA2C" (with to^ replaced by w 2 £,). The characteristic matrix, here denoted by 
Ch(H,to 2 ), is given by 


* The author has exercised the option to avoid printing all tabulated functions of 
frequency . 

~>;cC£D»WG PAGc BLANK NOT FILMED 


H 


Ch(H , to 2 ) 


" 1. O' 
0 1. 


where H, known as the dynamic matrix, is defined by 
H = M" 1 K , 


the matrices M and K being, respectively, the system mass matrix and stiffness 
matrix, that is [see equation (4) or (5)], 


M = 


m 0 


0 I 


, K = 


+ K 2 

K 2 L 2 

K 1 L 1 

- K 1 L 1 

k i l i 2 

+ K 2 L 


The equation formed by setting the determinant of Ch(H,ai ) to zero is a quadratic 
2 2 2 

in cj whose roots, and are the characteristic values (or eigenvalues of H) 

and also the squares of the undamped coupled natural frequencies. As the charac- 

teristic vector (or eigenvector) of H corresponding to u>. one may choose any non- 

1 2 

zero scalar multiple of either column of the adjoint of the matrix Ch(H,u>. _,), i = 1,2. 
Program TRROBM selects the second column of the adjoint and the reciprocal of the 
2,2 element as the scalar multiplier to get the vectors identified as "normalized" 
characteristic vectors on the output print. 


The vectors of the preceding paragraph could also be called modal columns of 
the "undamped" modal matrix. Notice that the modal columns have been "normalized" 
in a certain fashion, the "fashion" indicated. The author has not declared that such 
a normalization renders the modal matrix normalized with respect to the mass matrix. 
If it is desired that the modal matrix be normalized with respect to the mass matrix, 
one should select one of the four admissible expressions for the undamped modal 
matrix*, here denoted by 



’+11 

1 

T— H 


_+21 

^ 22 J 





*See Appendix A. 
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and post multiply it by the diagonal matrix 


n = 


n l 

0 


0 


n 


2 J 


where and a 2 are computed by 


"1 = (m *11 + 1 + 21 )_1/2 - "2 = <“ *u + 1 * 22 )_1/2 ' 

The expressions for n-^ and n 2 were found by simply demanding that the matrix 
V = $ n be such as to satisfy 


MASS = .10000000+01 MOMFNT OF INERTIA * .25000000*02 DELTA F = .50000000+ 00 
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BELOW WERE COMPUTED VJA THE TRAPEZOIDAL RULE APPLIED TO THE DEFINITE INTEGRALS DEFINING THEM 
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UNDAMPED COUPLED NATURAL FREQUENCIES : TCI . 999603GG+02 ( HZ ) f C2 - .20001981+03 (HZ) 
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o 

o 

CM 

O 

+ 

00 

CM 

o 

+ 

n 

CM 

o 

CT> 

C 

c 

T 

ro 

O 

+ 

0 

1 

CO 

CM 

O 

+ 

in 

CM 

O 

•f 

CM 

/•n 

CD 

CM 

o 

+ 

CD 

CM 

O 

+ 

<Ti 

CD 

co 

CM 

Oi 

co 


ro 


CO 

ID 

LO 

01 


CD 



O 

00 

CD 


CD 

CD 

CM 

r- 

in 

in 


" T— 

*- 

cd 

r- 

CD 

m 


CD 

ID 

r- 

, ~ 

CN 

0) 

in 

CD 

in 

O 

c 

00 

CD 


I s * 

CD 

f** 

r- 


C 

CD 

CM 

TT 

CD 


01 

CD 

in 



*- 

CM 

r- 

CD 

in 

CO 

ro 

l'. 


CM 

CM 


CO 

r- 

cm 


cm 

CO 

CM 

n 


ii 

ii 

n 

ii 

ii 

n 

ii 

i 


It 

n 

ii 

X 

XD 

XDD 

i- 

TD 

TDD 

- 

YD 

C 

>. 

YP 

Cj 

CL 

> 

YPDD 

li- 

u. 

L. 

Li_ 

LL 

LL 

LL 

Ll 

Ll 

Ll 

Ll 

ii 

o 

o 

O 

o 

c 

O 

c 

o 

O 

O 

O 

o 

LO 

l/l 

to 

to 

to 

to 

to 

LO 

to 

LO 

to 

to 

2 

s 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

a 

a 

a 

a 

Of 

Cl 

CL 

CL 

CL 

or 

a 

Of 



ro 

ro 

in 

o 

CD 

ro 

ro 

CO 

CM 

ro 


c 

G 

O 

O 

o 

O 

© 

O 

c 

O 

C 

o 

+ 

■+ 

+ 

■ 

+ 

4 

i 

+ 

•r 

1 


+ 

to 

00 

CO 

CD 

CD 

CM 

CD 

CD 

D* 

r- 


ro 

ao 

CM 

01 


O 

C 

CD 



CO 

TJ 


*r~ 

CD 

ro 

Ol 

CD 

in 


00 

CD 

CD 


ro 

00 

r- 

ro 

ro 

TT 

CD 

TT 

CO 

rr 

m 


CM 

C 



O 

in 

r- 

01 

CD 

'C 

Oi 

03 

CM 

■*— 

CD 


a- 

00 

CD 


CM 

in 


CD 

Tf 

CD 

in 


ro 

in 

•*- 

CD 

CD 


'LO 

n 

rj 

ro 

CD 

CD 

■*“ 

CD 

If.' 

D* 

CM 

CD 


in 




i» 

ii 

ii 

i. 

>' 

n 

If 

it 

Ii 

K 

II 

X 

C 


K 

O 

C 

> 

C 

O 

YP 

D 

0. 

c 

o 


X 

5 


K 

3 


**- 

G 


>- 

a. 



X 



1— 



> 




LL 

L- 

LL 

LL 

Ll 

Ll 

li- 

Ll 

li- 

LL 


LL 

o 

o 

O 

o 

O 

O 

CS 

O 

CS 

O 

c 

o 

L Li 

LL 

LL 

LJ 

LJ 

LJ 

LL 

LJ 

LU 

LL 

LL 

LJ 

Cl 

CL 

CL 

CL 

Or. 

CL 

CL 

a 

Of 

0L 

Of 

CL 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

<1 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

rs 

3 

O 

a 

O 

O 

O 

O 

a 

o 

a 

O 

o 

3 

to 

io 

CO 

CO 

c/; 

to 

to 

to 

to 

00 

01 

tn 

Z 

z 

-» 

Z 

Z 

Z 

z 

z 

z 

Z 

z 

z 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

LL 

LJ 

ll 

LL 

LJ 

LL 

LJ 

LJ 

LL 

jJ 


1 1 1 

2 

2 

gr 

2 

2 

2 

2 

2 


2 

£ 

2 


81 


UNDAMPED COUPLED NATURAL FREQUENCIES : FC1 = .93387255+02 (HZ) FC2 « .20317190+03 (HZ) 


* 







CD 


o 




■*“ 


CM 




< 


< 




C 


o 




UJ 


u 




2 


2 




Q 


c 




O 


c 




*— 


t— 




03 


CD 




z 


z 










Q 


o 




z 


z 




o 


o 




CL 


a. 




OO 


oo 




UJ 


LU 




a 


CL 




C£ 


CL 




o 


o 




CJ 


o 




CL 


CL 




o 


o 




K 


►— 




o 


o 


o 

CD 

UJ 


UJ 



CM 

> 


> 


< 

< 





o f" in 

O r-~ 

o 


U 


WOO 

wOO 

H 




2 + + 

2 + + 

1- 


t- 


O n ^ 

O n tt 

00 


00 


r- cm 

r- oo 

•— t • — ' 

oo 

1—1 i»— ‘ 

00 

0 — 00 

O - n 

CL OO 

Z 

CL 00 

Z 

h- Lf> 

K TT OC 

UJ UJ 

< 

UJ UJ 

< 

03 OC 

00 'C 

t- X 

MH 

*- z 


O ^ r 

O <T CO 

u u 

Q 

a u 

Q 

Z C N C 

Z CM CM 

< z 

< 

< z 

< 

— — in 

NH ^ 

cl w 

CL 

CL ~ 

CL 

o * • 

Q * • 

< W 


< w 


Z i • 

Z < 

X 


X 


o 

C 

V 


u 


Cl 

C- 

CM 

T- 


T— 

in r- in 

oo m in 

' O 

o 

' o 

o 

wOO 

wOO 

* + 

+ 

' + 

+ 

cl + + 

cl + + 

O ID 

c 

a in 

o 

cl — 

cl r- — 

uj 

o 

LU 'JD 

o 

O 03 OI 

O OUT- 

isi o 

o 

N <T> 

o 

u n id 

O hCJ 

oo 

o 

1-4 cn 

o 

00 ID 

in to 

-l CM 

o 

-J o 

o 

k ^ n 

t— 00 CO 

< r- 

o 

< •*- 

o 

2DO 

Z 't cn 

2 tt 

o 

2 O 

c 

H CN CD 

H O 

CL CM 


CL — 


o — n 

o in ^ 

O * 


C • 


“3 « ■ 

"3 • • 

z 


z 1 


O i « 

Q • 





< 

< 






82 






APPROVAL 


COMPONENT RESPONSE TO RANDOM VIBRATORY MOTION 
OF THE CARRIER VEHICLE 

By L. P. Tuell 


The information in this report has been reviewed for technical content. Review 
of any information concerning Department of Defense or nuclear energy activities or 
programs has been made by the MSFC Security Classification Officer. This report, 
in its entirety, has been determined to be unclassified. 


g. f. McDonough * 

Director, Structures and Dynamics Laboratory 


85 


TECHNICAL REPORT STANDARD II I Lb HAOe. 


1. REPORT NO. 2. GOVERNMENT ACCESSION NO. 

NASA TM -100309 


4. TITLE AND SUBTITLE 

Component Response to Random Vibratory Motion 
of the Carrier Vehicle 


7. AUTHOR(S) 

L. P. Tuell 


9. PERFORMING ORGANIZATION NAME ANO ADDRESS 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


12. SPONSORING AGENCY NAME ANO AOORESS 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


3. RECIPIENT’S CATALOG NO. 


5. REPORT DATE 

June 1987 


6. PERFORMING ORGANIZATION CODE 


8. PERFORMING ORGANIZATION REPORT W 


10. WORK UNIT NO. 

1 1. CONTRACT OR GRANT NO. 

13. TYPE OF REPORT & PERIOD COVERED 

Technical Memorandum 


14. SPONSORING AGENCY CODE I 


15. SUPPLEMENTARY NOTES 

Prepared by Structures and Dynamics Laboratory, Science and Engineering 
Directorate . 


16. ABSTRACT 

Two physical models of component plus supporting substructure are considered. 
Each model consists of a rigid body attached to a moving base by means of linear 
springs and viscous dampers. The second model differs from the first in that its 
dampers are elastically supported. The first model receives the more extensive 
treatment. Base motion, assumed a random translational motion parallel to a fixed 
axis, is prescribed only to the extent that the power spectral density (PSD) of its 
acceleration is given; and, as given, its plot on log-log graph paper is a series of 
straight line segments, each segment having an extremity in common with the 
adjacent segment. Closed expressions are given for the mean squares of base accel- 
eration, base velocity, and base displacement. The component is restricted to planar 
motion and allowed two degrees of freedom, one translational and one rotational. 
Integral expressions are given for the mean squares of component response variables, 
the transfer functions essential to mean square computation being available via the 
equations of motion. Closed expressions are given for mean squares of certain of the 
response variables for the case wherein the base acceleration PSD is constant. A 
very brief paragraph is given to stability of motion. 
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